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Abstract. We consider a charged particle moving in a static electromagnetic 
field described by the vector potential A(x) and the electrostatic potential V(x). 
We study the conditions on the structure of the integrals of motion of the first 
and second order in momenta, in particular how they are influenced by the gauge 
invariance of the problem. Next, we concentrate on the three possibilities for inte- 
grability arising from the first order integrals corresponding to three nonequivalent 
subalgebras of the Euclidean algebra, namely (Pi,P 2 ), (Z/ 3 ,P 3 ) and (Li, L 2 , L 3 ). 
For these cases we look for additional independent integrals of first or second or¬ 
der in the momenta. These would make the system superintegrable (minimally 
or maximally). We study their quantum spectra and classical equations of mo¬ 
tion. In some cases nonpolynomial integrals of motion occur and ensure maximal 
superintegrability. 

Keywords : integrability, superintegrability, classical and quantum mechanics, 
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1. Introduction 

The purpose of this article is to initiate a systematic study of integrable and 
superintegrable systems in the presence of a magnetic field in three-dimensional 
Euclidean space E 3 , both in classical and quantum mechanics. 

We recall that a classical Hamiltonian system with n degrees of freedom is 
(Liouville) integrable if it allows n integrals of motion X \,..., X n including the 
Hamiltonian. They must be in involution, be well defined functions on phase space 
and be functionally independent. The system is superintegrable if further integrals 
Y a exist with 1 < a < n — 1. They must also be well-defined functions on phase 
space, the set { Xi,Y a } must be functionally independent, however the integrals Y a 
need not be in involution with X,; nor with each other. On the contrary, the set 
of all integrals generates a non-Abelian algebra under Poisson commutation. In 
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quantum mechanics the definitions must be modified. The integrals X\,... ,X n will 
be well-defined operators in the enveloping algebra of the Heisenberg Lie algebra, i.e. 
polynomials in the coordinates x 3 and pi , or more general objects, such as convergent 
series in these objects. Instead of functional independence we shall require polynomial 
independence. In other words no nontrivial fully symmetrized polynomial (Jordan 
polynomial) in the integrals of motion should vanish. 

The best known superintegrable systems in E n (n > 2) are the Coulomb-Kepler 
system and the isotropic harmonic oscillator with their o(n + 1 ) and su(n) symmetry 
algebras, respectively [D13 Gam. It follows from Bertrand’s theorem [5] that these are 
the only spherically symmetric maximally superintegrable systems in E 3 (and actually 
in E n (n > 2)). Most of the recent research on superintegrability concentrated on 
“natural” Hamiltonians with scalar potentials. For a recent review see |!|. For early 
systematic work on superintegrability in E 2 see [TUB] and in E 3 see also 

mmm and references therein. 

Earlier work on integrability with magnetic fields (velocity dependent potentials) 
mainly concerned the two-dimensional case with notable 

exceptions [22] [23] . 

A particle moving under the influence of a static electromagnetic field is described 
by the classical Hamiltonian 

H =^{p + A(x)) 2 + V(x) (1) 

where p = (pi,P 2 ,P 3 ) are components of the linear momentum, x = (x±, x 2 , X 3 ) = 
( x,y,z ) the spatial coordinates, the vector and electrostatic potential A(x), V(x) are 
functions only of the coordinates x and we choose the units in which the mass of the 
particle has the numerical value 1 and the charge of the particle is —1 (having an 
electron in mind as the prime example). 

We recall that the equations of motion of the Hamiltonian Q are gauge invariant, 
i.e. that they are the same for the potentials 

A'(x) = A(x) + V X , V'(x) = V(x) (2) 

for any choice of the function x(A)- Thus, the physically relevant quantity is the 
magnetic field 

- - dAi 

B = V x A, i.e. B j = ej k i-^~ (3) 

rather than the vector potential A(x) (where ejki is the completely antisymmetric 
tensor with £123 = 1 ). 

We shall also consider the quantum Hamiltonian defined as the properly 
symmetrized analogue of |lj) in terms of the operators of the linear momenta P 3 = 
—ihd/dxj and coordinates Xj — x 3 : 

h = \Y. + p M*) + M*) p i+Mi) + M)- ( 4 ) 

j 

The operators Aj(x) and V(x) act on wavefunctions as multiplication by the functions 
Aj(x) and V(x), respectively. 

On the quantum level, the gauge transformation ([ 2 ]) demonstrates itself as a 
unitary transformation of the Hilbert space. Namely, let us take 

Utp(x) = exp f tX(x) ) ■ V>(z)- 


(5) 
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Applying ([5]) on the states and the observables we get an equivalent description of the 
same physical reality in terms of 

%j) = Uijj, 6^6' = UOUK (6) 

In particular, the following observables transform covariantly 

(Pj + As) -+ u(Pj + Aj)u t = Pi + Aj, v -+ uvtf = v. (7) 

We recall that the quantum dynamics may not depend only on the magnetic held B. 
In nontrivial topologies of the configuration space, e.g. when singularities are present 
in the potential, not all vector potentials A inducing the same magnetic held B are 
gauge equivalent, as the observable Aharonov-Bohm effect demonstrates [24] . 


2. The conditions for the integrals of motion 


Let us consider integrals of motion which are at most second order in the momenta. 
Because our system is gauge invariant ©> 0 we hnd it convenient to express the 
integrals in terms of gauge covariant expressions 

Pj = Pj + A j > Pf = Pj + Aj (8) 

rather than the momenta themselves. The operators ([8]) do not commute amongst 
each other. They satisfy 

[Pf, Pk ] = -itejkiBu [P/, X k ] = —iHSjh (9) 

and analogous relations for Poisson brackets in classical mechanics. 

Classically, we write a general second order integral of motion as 

3 3 1 3 

x = '52 h j(x)pfpf + H ^\tjki\nj(x)p£pf + J2sj(x)pf + m(x). (10) 

j= 1 j,k,l=l J=1 


We dehne the Poisson bracket in the standard manner 


{a(x,p),b(x,p)}p. B . = 


j =i 


da db _ db da 
dxj dpj dxj dpj 


( 11 ) 


The condition that the Poisson bracket of the integral (10) with the Hamiltonian (|T]) 
vanishes 


{H,X} p . b .= 0 (12) 

leads to terms of order 3, 2,1 and 0 in the momenta and respectively to the following 
equations: 

Third order 


d x h\ = 0, dyhi = -d x n 3 , d z hi = -d x n 2 , 

^xb‘2 dyTiQi dyh 2 0, d z h 2 dyTi\- ) 

d x h$ — d z n 2 , dyh$ d z n\^ d z h$ — 0, 
V • n = 0. 


(13) 
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d x s i = n 2 B 2 - n 3 B 3: 
d y s 2 = n 3 B 3 - n 1 B 1 , 
d z s 3 = n\B\ - n 2 B 2 , 

dySi + d x s 2 = n\B 2 - n 2 Bi + 2{hi - h 2 )B 3 , (14) 

d z S\ + d x s 3 = n 3 Bi — n\B 3 + 2 (h. 3 — hi)B 2 , 

<9 y S3 + = ri 2 l ?3 — n 3 B 2 + 2(h 2 — h 3 )B\. 

It follows that we have 

V • s = 0. 


First order 

d x m = 2hid x V + n 3 d y V + n 2 d z V + s 3 B 2 - s 2 B 3 , 

d y m = n 3 d x V + 2h 2 d y V + + s 3 B 3 - s 3 B ± , (15) 

d z m = n 2 d x V + nid y V + 2h 3 d z V + s 2 Bi - siB 2 . 

Zero order 


s-W = 0 . 


(16) 


Equations (13) are the same as for the system with vanishing magnetic field and 


their explicit solution is known. They imply that the highest order terms in the 
integral ([To]) are linear combinations of products of the generators of the Euclidean 
group pi , p 2l p 3 i li, l 2: l 3 where lj = i^jkiXkPi (up to a redefinition of lower order 
terms). When expressed explicitly in terms of the covariant expressions (18]), we have 


X= a abYa Y b A + S i ( f ) Pf + m (®) 

1 <a<b<6 j =1 


where Y A = {pi,p 2 ,p 3 ,li,l 2 ,l. 


3 )i 


l t = J2 


€jkl%kPl 


(17) 


(18) 


and a a b £ R- The functions h, n are expressed in terms of the constants a a b as follows 
hi = a 66 y 2 + (~a 56 z - aie)y + a 55 z 2 + a 15 z + an, 

h 2 = a e ex 2 + (— a 4 6 z + a 26 )x + 044 z 2 - a 2 4 Z + a 2 2 , 

h 3 = a 55 x 2 + (-a 4 5 y - 0:35)2; + 044 y 2 + o 34 y + o 33 , ( 19 ) 

tii = ~a 56 x 2 + (a 46 y + 0452: - a 2 5 + o 36 )a; + (-2a 44 2 + a 2i )y - a 34 2 + a 2 3, 
n 2 = (a 56 y - 2a 55 2 - 045)2; - a 4 6 y 2 + (0452 + o 36 + a 14 )y + 0352 + ai 3 , 
n 3 = (~2a 66 y + ai 6 + o 56 2)2; + (a 46 2 - a 26 )y - a 45 2 2 + (a 2 5 - 044)2 + a 12 . 

Since we have the relation 

p-T= 0 ( 20 ) 


one of the parameters 044 , 025,036 is irrelevant and can be set equal to any chosen 
value. However that explicitly breaks the symmetry of our equations; thus we prefer 
to keep all parameters 044 , 025,036 in our formulas and we use ( 20 ) as a constraint. 
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The conditions (14) allow us to express all the second order derivatives of s 
explicitly in terms of ft, h, B and their derivatives 

9 2 x s\ = 9 x n 2 B 2 + n 2 d x B 2 — d x n 3 B 3 — n 3 d x B 3 , 

9 xy s i = d y n 2 B 2 + n 2 d y B 2 — d y n 3 B 3 — n 3 d y B 3 , 
d xz si = d z n 2 B 2 + n 2 d z B 2 — d z n 3 B 3 — n 3 d z B 3 , 
dy y si = n\d x Bi - 3d x n 3 B 3 - n 3 d x B 3 + d x n A B { + 2h 1 d y B 3 - d y n 2 Bi 
- n 2 d y B 1 + n\d y B 2 - 2h 2 d y B 3 + d y n 3 B 2 , 


3 13 13 

9 yz si = — - d z n 2 Bi — -n 2 d z B\ + -d y n 3 Bi + -n 3 d y Bi — -d x n 2 B 3 

- \n 2 d x B 3 + h\d z B 3 + ]^d y n 1 B 3 - h 2 d z B 3 - ]^d z niB 2 + ^ nid z B 2 

13 1 

- -n\d y B 3 + -d x n 3 B- 2 - hid y B 2 + h 3 d y B 2 + -n 3 d x B 2 - h 2 d x B x + h 3 d x B 3 , 


d 2 zz si = 3d x n 2 B 2 + n 2 d x B 2 - d x niB 1 - nid x B 3 - 2hid z B 2 - d z n 3 B 3 

- n\d z B 3 + n 3 d z B 1 + 2h 3 d z B 2 + d z n 3 B x , 

9l x S2 = - d y n 2 B 2 - n 2 dyB 2 + 3d y n 3 B 3 + n 3 d y B 3 + 2h 1 d x B 3 - d x n 2 B 1 

- n 2 d x B 1 + Ul d x B 2 — 2h 2 d x B 3 + d x n\B 2 , 
dl y s 2 = d x n 3 B 3 + n 3 d x B 3 - d x n 3 Bi - n 3 d x B 3 , 

9l z S2 = ^d z n 2 B 1 - ^ n 2 d z Bi - ^ d y n 3 B 1 - ^ n 3 d y B 1 - ^d x n 2 B 3 + ^n 2 d x B 3 

3 3 11 

+ h\d z B 3 + -d y n\B 3 - h 2 d z B 3 + -d z mB 2 + -nid z B 2 + -nid v B 3 

3 1 

- ~j9 x n 3 B 2 + hid y B 2 - h 3 d y B 2 - -n 3 d x B 2 + h 2 d x Bi - h 3 d x B lt 


9yyS2 9yTl 3 B 3 ~\~ Tl 3 dyB 3 0y7l\B\ Tl\dyB\, 

d yz s 2 = d z n 3 B 3 + n 3 d z B 3 - d^iBx - n 1 d z B 1 , ( 21 ) 

d zz s 2 = — nidyBi + d y n 2 B 2 + n 2 d y B 2 - 3d y niBi + 2h 2 d z Bi - d z n 3 B 2 

- n 3 d z B 2 - 2h 3 d z B 1 + d z n 2 B 3 + n 2 d z B 3 , 

d xx s 3 = n 3 d z B 3 - 3d z n 2 B 2 - n 2 d z B 2 + d z n 3 B 3 - 2hid x B 2 - d x mB 3 

- n\d x B 3 + n 3 d x Bi + 2h 3 d x B 2 + d x n 3 Bi, 

3 1113 1 

9l y s 3 = -d z n 2 Bi + -n 2 d z Bi - -d y n 3 B 1 + -n 3 d y Bi + -d x n 2 B 3 + -n 2 d x B 3 

3 3 11 

- h\d z B 3 - -d y n\B 3 + h 2 d z B 3 - -d z mB 2 - -nid z B 2 - -nid y B 3 

+ ^d x n 3 B 2 - h\d y B 2 + h 3 d y B 2 - ^n 3 d x B 2 + h 2 d x Bi - h 3 d x B 1 , 


9 xz s 3 = — d x n 2 B 2 — n 2 d x B 2 + d x n\B\ + n\d x B \, 

9ly s 3 = - d z n 3 B 3 - n 3 d z B 3 + 3d z n 3 Bi + nid z Bi + 2h 2 d y B 1 - d v n 3 B 2 
— n 3 d y B 2 — 2h 3 d y Bi + d y n 2 B 3 + n 2 d y B 3 , 

9 yz s 3 = — d y n 2 B 2 — n 2 d y B 2 + d y n\Bi + n\d y Bi, 

9 2 zz s 3 = - d z n 2 B 2 - n 2 d z B 2 + d z niBi + nid z Bi. 
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Taking various first and second order derivatives of these expressions and comparing 
them, e.g. 

d v( d xx s i) = dx{d 2 xySj ), 

we arrive at compatibility conditions for h, n and B. Altogether there are 6 
independent second order PDEs relating them which can be solved e.g. for 

dyyBi , d zz Bi, d xx B 2 , d zz B2 , d xx B 3 , d yy B 3 

but they are too cumbersome to reproduce here. 


In the quantum case we have to consider a properly symmetrized analogue of (10). 
We choose the following convention 

1 = p/p/} + + +m ^ ( 22 ) 

j=i j,k,l=i j=i 

where { , } denotes the symmetrization 


{a, b} = - (ab + ba'j 


and hj(x),nj(x), Sj(x),m(x) are to be interpreted as the corresponding operators 
of multiplication by the given function. All possible choices of symmetrization are 
equivalent up to redefinition of the lower order terms. 


Imposing the condition that the integral of motion in the form (22 1 commutes 
with the Hamiltonian Q we obtain a similar set of conditions as above. In particular, 
the conditions (131 remain the same, i.e. their solution has the same form (191. The 


conditions (14) get apparent quantum corrections 
d x si = n 2 B 2 - n 3 B 3 , 
d y s 2 = n 3 B 3 - n\B\ : 
d z s 3 = n\B\ - n 2 B 2 , 
d y s 1 + d x s 2 = n 3 B 2 — n 2 B 3 + 2 (hi — h 2 )B 3 + ihd 
dzSi + d x s 3 = n 3 B x - n 3 B 3 + 2 (h 3 - h 3 )B 2 + ihd. 


xx n 3> 
|2 


d v s 3 + d z s 2 = n 2 B 3 - n 3 B 2 + 2 (h 2 - h 3 )B 3 + ihd y n 3 . 


Due to the explicit solution (19) we have 


(23) 


^xx^3 0 ; 


— 0, 


d yy n 1 = °> 


i.e. these apparent quantum corrections vanish and equations (231 and (14) are the 
same. 

Similarly, in the quantum version of equations (15) we obtain some extra terms 
but they vanish once the solution (19) of equations (131 is substituted in, i.e. also the 
conditions (15) remain in the quantum case. 

The situation is however different for the last equation (16) which does obtain 
fi 2 -proportional corrections. Using equations (13), (14), (15) and their differential 
consequences we find 


s- VU + — ( d z nid z Bi - d y n 1 d y B 1 + d x n 2 d x B 2 - d z n 2 d z B 2 + (24) 
+d y n 3 d y B 3 - d x n 3 d x B 3 + d x n 3 d y B 2 - d y n 2 d x Bi) = 0. 
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Notice that the last line seems to violate the obvious symmetry of our problem under 
Euclidean transformations. However, this is not the case in view of the identity 

d x n\dyB 2 - d y n 2 d x Bi = d y n 2 d z B 3 - d z n 3 d y B 2 = d z n 3 d x B\ - d x n 1 d z B 3 (25) 

which is a consequence of 

V • B = 0, V • n = 0. 


Notice that for the special case of first order integrals the conditions (13 H16) 
become significantly simpler. Namely, the conditions (131 do not arise. The right hand 
sides of the conditions (141 vanish and thus equations 13 imply that the first order 
term Xq=i s j(.%)pf i n the integral is a constant linear combination of the covariant 
linear and angular momenta p 3 ,p 2 ,p 3 , l 3 , l 2 , l 3 (see 0, ©)■ The conditions ( fl5| ) 
simplify to 


d x m = s 3 B 2 - s 2 B 3 , d y m = s 3 B 3 - s 3 B 1} d z m = s 2 B ± - s 3 B 2 (26) 

and imply first order compatibility conditions relating B and s 


d y s 3 B 2 + s 3 d y B 2 - d v s 2 B 3 - s 2 d y B 3 - d x s 1 B 3 - s 1 d x B 3 + d x s 3 B 1 + s 3 d x B 1 = 0, 
d z siB 3 + s\d z B 3 - d z s 3 Bi - s 3 d z B 3 - d y s 2 Bi - s 2 d v B x + d y s 3 B 2 + s 1 d y B 2 = 0(27) 
9 z s 3 B 2 + s 3 d z B 2 — d z s 2 B 3 — s 2 d z B 3 — d x s 2 Bi — s 2 d x Bi + d x S\B 2 + s\d x B 2 = 0. 


The condition (161 remains the same as for the second order integral. In this case it 
gets no quantum correction, i.e. (16) and (24) now coincide. 


Let us now turn our attention to the situation where the Hamiltonian 0 or § 
is integrable in the Liouville sense, with at most quadratic integrals. That means 
that in addition to the Hamiltonian itself there must be at least two independent 
integrals of motion of the form (10) or (22) which commute in the sense of the 
Poisson bracket or Lie commutator, respectively. Independence is to be understood as 
functional independence in the classical situation and in the sense that no nontrivial 
fully symmetrized polynomial in the given operators vanishes in the quantum case. 

Since the highest order conditions (13) are the same whether or not there is 
a magnetic field present, the first step of the analysis can be performed as in [5j, 
leading to 11 nonequivalent possibilities for the functions h,n. Next one should 
look into equations (141, (151, ( [l6| (resp. ( [24| ) and their consequences to determine 
the nonequivalent possible choices of the magnetic field B(x) and the electrostatic 
potential V{x). This is the approach used by A. Zhalij in [22] for the special case of 
the integrals X\ = Pf + ..., X 2 = P| + ... 


We shall follow a different route here. Keeping in mind that our main goal is 
to arrive at examples of superintegrable systems with nonvanishing magnetic held 
we shall assume that the integrability arises in the simplest way possible. Namely, 
we assume that there are at least two independent first order integrals for our 
Hamiltonian. 

Assuming that we have 


x i = i\^i + PiPi + m^x), X 2 = + fi l 2 pf + m 2 {x) (28) 


(or its quantum analogue) we may use the Euclidean transformations to simplify 
A'i, X 2 . Another allowed transformation is replacing X \, X 2 by an arbitrary regular 
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Xi —y X\ — K\X\-\- K2X2, 


X2 —y X2 — Ai-X"i + X2X2, 


det 


«i 

Ai 


K 2 

A 2 


^ 0 . 


For convenience, we redefine the yet unknown functions toi(:e), m 2 (af) as needed 
without renaming them. 

We arrive at the following possibilities 

• If we have 71 = 72 = 0 then we can set X\ and X 2 by rotation and linear 
combination to 

X\ = Pi + TOr(x), X 2 =pi + m 2 {x). (29) 

• If (71,72) 7 ^ (0,0) we can transform e.g. X\ by rotation and translation into 
Xl = if + /3pf + nii (x) ■ 

— Assuming that the integrability arises directly at the first order, i.e. that 
{Ad,X 2 } p.b. = 0 , we arrive at a single possibility 

X 1 =lf + na(x), X 2 =ps +m 2 (x). (30) 

— However, there is another option - to allow X\ and X 2 to be not in involution 
and expect the second commuting integral to arise via Poisson brackets and 
polynomial combinations of Xi,X 2 . Thus we may up to rotation and linear 
combination take 

Xi = if + fipf +mi(x), X 2 = ali + fi l 2 pf + m 2 (x), er = 0,1.(31) 
In order to have nontrivial dynamics, i.e. nontrivial electric and/or magnetic 
field, we cannot have the full Euclidean algebra represented in terms of the 
integrals of motion. Thus we must require that the algebra generated by 
the highest order terms I 3 + ftps and al 1 + fi l 2 pi in (311 via Poisson brackets 
closes as a proper subalgebra of the Euclidean algebra. We have the following 
options: 

(i) The algebra isomorphic to su(2) 


1 = t 3 - 1 - rm(x), X 2 = if + m 2 (x), 


A-, = fi 

X 3 = {Xi,X 2 }p.b. = l 2 +m3(x). 
(ii) The algebra isomorphic to the Euclidean algebra c 2 

M '-pi 


(32) 


Xi = I 3 +P 3 +mi{x), 


X 2 = pf 


m 2 (f), 

X 3 = {X -[, X 2 }p.b. = p 2 + rri 3 (x). 

This case is, however, already included in (29) as a special subcase. 


3. Superintegrability for the integrable system with integrals Pi,P 2 


Let us start our detailed investigation by considering the case of the integrals (29) 
Xi=pi+mi(x), X 2 =pi + m 2 (x). 


The condition that Xi and A' 2 are in involution is equivalent to 


d v mi - d x m 2 = B 3 . 


Equations (15) reduce to 


d x mi = 0, dyUii = B 3 , 
d x m 2 = -B 3 , d v m 2 = 0, 


d z mi = -B 2 , 
d z m 2 = Bi 


(33) 


(34) 
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and together with (331 imply that 

B 3 {x) = 0 , Bj(x) = Bj(z), mj(x)=mj(z), j = 1,2. 

Writing the components of the magnetic field as 

B j (x)=F'(z), J = 1,2, (35) 

we have 

m j(x)= F j( z )’ j = 1 ,2. (36) 

We choose a suitable vector potential in the form (satisfying the Coulomb gauge 
condition VA = 0) 

A(x) = (F 2 (z),-F 1 (z), 0) (37) 

and from the conditions ( fl6| ) we find that 

V(x) = V(z). (38) 

Plugging all the information obtained about functions A, B, mj into the assumed form 


of the integrals (29) we find a very simple solution (unique up to the choice of gauge) 
Xi = pi, X 2 = p 2 . (39) 

The same result arises also in the quantum case, via essentially the same arguments. 

Let us now assume that our system ([Tj) with the potentials 

A{x) = (F 2 (z),-F 1 (z),0) , V(x)=V(z) 

is superintegrable, i.e. that an additional independent integral of motion exists. For 
simplicity, let us assume that it is of first order in momenta. Up to addition of Xi 
and X 2 we have 

X 3 =flf+ppi + m 3 (x). (40) 


We consider equations (151 and their compatibility conditions which take the form 
12XF[' - 7l yF" - p F” - 73 F' = 0, 7l U' - j 2 F[ = 0, 

- 72 ^ , + 7 i 2 /F"- 73 F(+^F" = 0. (41) 

We arrive at two distinct possibilities: 

• If 7 i + 7 f 7^ 0 then F" = F." = 0, i.e. the magnetic field (35) is constant. This 
case has already been well studied in the literature, see e.g. [23 ESI- Solving 


equations (15) and (16) we find that the electrostatic potential is constant too, 


i.e. we have a motion in constant magnetic field and no electric field. This system 
is superintegrable and exactly solvable as follows. Without loss of generality we 
can rotate the coordinate system so that we have 

B(x) = (B,0,0), A(x) = (0,-Bz,0), V(x) = 0. (42) 

Four independent first order integrals exist in this case. We write them down in 
the classical situation, quantum mechanically they are the same expressions in 
terms of operators. They read 


X\ — Pi , 


X 2 — p 2 , 


-Y 3 =p 3 - By, 


x 4 — ix + — y 2 ^- 


(43) 
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The Hamiltonian can be expressed in terms of X\,, X 4 as 
H = ~{p\ + p\ + p\) - Bzp 2 + — z 2 = - (Xf + X^ + X%) + BX 4 . ( 44 ) 

The classical equations of motion 

X=P 1 , y=P 2 -Bz, Z=p 3 , ( 45 ) 

Pi =0, p 2 = 0, p 3 = B(p 2 - Bz) 

are solved explicitly as 

X{t) = XQ+Pit , 

y(t) =yo-j^+ cos (Bt)^ - sin (Bt) (z 0 - P ^j , (46) 

z (t) = P g + + cos (Bt) (^zo - , 

Pi(t) = Pi, P 2 {t) = p°, p 3 (t) = cos {Bt)p% + sin {Bt){pl - Bz 0 ) 


where ( Xo,yo,Zo ) are the initial coordinates and (Pi,P 2 ,P 3 ) the initial momenta. 
From (46) we find that the trajectory is a helix with axis parallel to the a;-axis and 
the integrals X 3 , X 3 , X 4 determine the diameter and position of the enveloping 
cylinder in the yz- plane 



(47) 


When pi = 0 the helix collapses into a circle in the plane x = xq. Thus the 
problem reduces to the two-dimensional one. Let us now restrict to the case 
Pi 7 ^ 0 . 

By inspection of the solution of the equations of motion one finds that this system 
is maximally superintegrable with, however, the fifth independent integral not 
polynomial in momenta. It reads 

X 5 = (Bz - p 2 ) cos - P 3 sin ■ ( 48 ) 

How to interpret this integral in the quantum case is not too clear. However, if we 
restrict ourselves to the subspace in the Hilbert space defined by the constraint 

Pifi(x) = ki%j}(x), fci/0 (49) 

we can expand X 5 into a convergent Taylor series in . We can then interpret 
the quantum integral of motion X$ as an operator in the “extended enveloping 
algebra” of the Heisenberg algebra. 

Alternatively, for the classical Hamiltonian we can perform a canonical 
transformation in the x,pi plane 
2 

Pi = x=—, Pi>0 (50) 

2 pi 

which transforms the integrals into 

H = Pi + \((p 2 - Bz) 2 +p\) , Xi=pi, X 2 =p 2 , 

X 3 = p :i - By , X 4 = h + j (z 2 - y 2 ) , 

X 5 = (Bz — P 2 ) cos (Bx) — p 3 sin ( Bx ). 


(51) 
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Thus X 5 becomes a first order polynomial in the momenta. The price is that H 


no longer has the “natural” form (44). 


The integrals X 4 ,..., X 5 give rise to a Lie algebra of integrals of motion in the 
following manner. We define additional functionally dependent integrals 

X 6 = {X 4 , X 5 }p,B. = (j >2 - Bz) sin - P 3 cos > X 7 = 1 (52) 

and redefine the first integral to be 

If 2 n 2 

X ' = -f = I' ^ 

The Poisson brackets now form a 7-dimensional Lie subalgebra of integrals of 
motion in the algebra of observables on our system, with the Lie brackets 


{ , }p.B. 

Xi 

X 2 

X 3 

X 4 

X 5 

X 6 

X 7 

Xi 

0 

0 

0 

0 

-BXe 

bx 5 

0 

X 2 

0 

0 

bx 7 

-X 3 

0 

0 

0 

X 3 

0 

-bx 7 

0 

X 2 

0 

0 

0 

X 4 

0 

X 3 

-X 2 

0 

Xq 

X 5 

0 

X 5 

bx 6 

0 

0 

-X 6 

0 

-bx 7 

0 

X 6 

-bx 5 

0 

0 

X 5 

bx 7 

0 

0 

x 7 

0 

0 

0 

0 

0 

0 

0 


(54) 


This algebra is solvable with 5-dimensional nilradical spanned by X 2 , X 3 , X 5 , Xq, X 7 . 
The nilradical is isomorphic to the Heisenberg algebra in two spatial dimensions 
(n 5; 3 in the notation of [27]). The element X 7 spans the center. Its Casimir invari¬ 
ants are the central element A ' 7 and two second order invariants 2X| X 7 +X|+Xfi, 
2(13X4 + Xi)X 7 + Xf + X| which both reduce to the Hamiltonian (44) once the 
explicit form of the integrals is inserted into them. 

We observe that in terms of the original integral X- t instead of its square X 7 we 
obtain an infinite-dimensional loop algebra. 

As it is well-known in the physics literature (see e.g. [25] . pg. 220), the Schrodinger 


equation for the Hamiltonian (44) separates in Cartesian coordinates as follows 


ip(x) = f(z) exp ( -kix 


exp ( - k 2 y 


(55) 


h 2 f(z) = ((Bz - k 2 ) 2 +k\- 2 E) f(z), 

Xi ip(x) = k\ip(x), X 2 ip(x) = k 2 ip(x). 

The reduced Schrodinger equation ( |55| ) is the stationary Schrodinger equation 

fc 2 

for the 1-dimensional harmonic oscillator with the energy E — frequency 
oj = B and the center of the force at z = k 2 /B. Thus the spectrum of the 


Hamiltonian (44) is continuous due to the arbitrary momentum k\ and reads 

1 


E = 


k 2 


hB 


n e No, fci £ 


2 V 2 / 

The eigenvectors are expressed in terms of Hermite polynomials 


(56) 


^n,ki ,&2 (*^) — H-i 



exp 


- (kxx + k 2 y) 



eXP| ^2 h 
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It was conjectured in [28] that all maximally superintegrable quantum systems in 
Euclidean spaces are exactly solvable. Equations (56 1 , (57) confirm this conjecture 
for a particle in a constant magnetic field. This is true even though one of the 
independent integrals X 3 is not polynomial in its natural setting. 

• If 71 = 72 = 0 we have to distinguish further subcases. If 73 = 0 the system 
again reduces to a constant magnetic field and vanishing electric field. If 73 ^ 0 
and /? = 0 the magnetic field must vanish and we are in a situation without 
vector potential, which is not of interest here (and already well studied, see e.g. 
El Ml (III S]). However, if 73 7 ^ 0 , j3 / 0 (and without loss of generality we can 
assume 73 = 1) we obtain a nontrivial solution for A and V 


A(x) = ( —A cos 


</>o 


B (x) = cos 
V(x) = const. 


, — A sin 


+ 0o \ A . / z + 0o 


■ 00 


,0 


’ ~p Sm 


P 


,0 


(58) 


where A > 0 and 0o are integration constants. We can simplify it by Euclidean 
transformations and a shift of the potential to 


A(x) = ( — A cos ( ^ ) , — A sin ( A ) ,0 


0 


P 


B{x) = „ cos ( — , — sin — , 0 


A 


V(x) = 0. (59) 


PJ’P \P 

The integral of motion X 3 (40) reduces to 

X 3 =l 3 + fip 3 (60) 

in the gauge chosen above. As before, this calculation is the same in both classical 
and quantum mechanics. 

The classical equations of motion take the form 


A 


P 


Pi(t)=0, p 2 (t)=0, p 3 (f) = — -sin —- pi(f)+cos —- p 2 (i) 


z(t) 


P 


z{t) 


x(t) = piit) — A cos 


z{t) 


y(t) =P 2 (t) - A sin 


z(fi 0 
P 


P 


z{t)=p 3 (t) (61) 


and can be solved by quadratures. Namely, we express the conserved momenta 
in the polar form 

(62) 


Pl =p cos \ — 


P 2 =P sm [ — 


where p > 0 and 0 P are constants, and find a second order equation for z(t) 


Ap . 

zip ) = — y sm 


z{t) - 0 r 


The order of this equation can be lowered, obtaining 


\ ( z(t)f = Ap ( cos 


Z{t) - 0 f 


K > — 1 . 


(63) 


(64) 


(k < —1 is unphysical since then (64) doesn’t have real solutions.) We substitute 
zip) = 0 P + /? arccos (C(^)) (65) 
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and we find a separable first order ODE for ((t) 

mr = - ixcw + ixcw+«) (66) 

solvable in terms of elliptic integrals. Explicitly, we change the independent 
variable 


t = . A T 

y/2Ap 

to have a simpler equation 

(C(r)) 2 = -(C(r)^l)(C(r) + l)(C(r) + K ). 


(67) 

( 68 ) 


The solution depends on the value of the integration constant k, namely 1 < n < 
— 1 and 1 < k demonstrate different behavior, and on the initial value for CM- 
For K>lwe find a solution in the form 


CM 


_ l-K 2 _ 

2 sn 2 (iyX + T(r - r 0 ), -K-l 


(69) 


For — 1 < k < 1 we find a solution in the form 

2(1 - k) 


CM = 


2 — (k + l)sn 2 - r 0 ), 


- 1 . 


(70) 


The equations of motion for x(t),y(t) now reduce to quadratures (611 in terms 
of the Jacobi elliptic function sn. Solving them numerically we obtain the 
trajectories for our system. For — 1 < k < 1 they are bounded in the plane 
perpendicular to ( p \, P 2 ,0) and appear like a deformed helix whose axis is parallel 
to the vector (pi,P 2 > 0). For 1 < n they are no longer bounded in the ^-direction 



Figure 1. Sample trajectory for —1 < k < 1 (with A = 3, 0 = 3,pi = 1,P2 = 
0, z(0) = 0.08, y{ 0) = 0.05, z(0) = 0, i(0) = 3.2) 


and appear like a deformed helix whose axis is no longer parallel to the ccy-plane. 
The value n = 1 appears to be a limiting case of the n > 1 situation. 

In the quantum case the stationary Schrodinger equation 

Hfi(x) = Efi(x) 


separates in Cartesian coordinates. We have 


K 


'C’M = xM ex P i“F cos (^ff) x ex P i“T sin (^fc) y 


I< 


Xii/j(x) = K cos (fiM'C’M, 


X 2 ip{x) = K sm((j> K )ip(x), 


h 2 xM = ( -2AK cos ( A - (j) K ) + A 2 + K 2 - 2 E ) xM- 


(71) 
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Figure 2. Sample trajectory for re > 1 (with A = 3, 0 = 3, pi = 1,P2 = 0,x(0) = 
0.08, y(0) = 0.05,2(0) = 0.1, i(0) = 3.5) 



Figure 3. Sample trajectory for re = 1 (with A = 3, 0 = 3, pi = 1,P2 = 0,x(0) 
' y(0) = 0.05, 2(0) = 0, 2(0) = 2y/3) 


The separated equation for \{ z ) is solved in terms of Mathieu sine and cosine 
functions C(a,q,x) and S(a,q,x ), i.e. independent solutions of the Mathieu’s 
differential equation 

y + {a — 2 qcos( 2 x))y(x) = 0 , 


as follows 


<t>k 


X (z) = C\ C ( -4^ (A 2 + K 2 - 2E), -A^AK, ^ - ^ ) + 


H 2 ' 


2(3 


(72) 




This leads to the conjecture that the Hamiltonian system with the potentials (59) 
is maximally superintegrable but an explicit calculation shows that that 
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hypothetical further integral cannot be of order one or two in momenta. In 
the classical case an additional fifth integral can be found from the Hamiltonian 
flow ( f6T| ) . We use the method of characteristics to arrive at the equation 
d y dz 


P '2 A sin (-3 


(h) ^2Acos(^)pi + 2 Hsin p 2 + ■ 


(73) 


where u = p 2 — 2A (pi cos(|)) + p 2 sin f = 2H — Xf — X 2 — A 2 is a constant 
of motion. Its solution is expressed in terms of Jacobi elliptic functions whose 
arguments depend on the momenta pi,p 2 and p 3 . Hence is not polynomial in the 
momenta at all. 


4. Superintegrability for the integrable system with integrals L 3 ,P 3 


Let us perform a similar analysis for the case (30) 

X 1 =l 3 +m 1 (x), X 2 =pf+m 2 (x). 

Requiring that X- t and X 2 are in involution we find the condition 
xBi + yB 2 + xd y m 2 — yd x m 2 — d z m\ = 0. 


Equations (15) reduce to 

d x rn-[ = -xB 3 
d x m 2 = B 2 , 


dyirii = -yB 3 , cLtoi = yB 2 + xB u 

d y m 2 = -Bi, d z m 2 = 0. 

Solving these equations and their compatibility conditions we find that 

m\{x) = —F 2 (R), m 2 (x) = Fi(R), R= ^x 2 + y 2 , 


Ax) = (-±F 2 {R), ^F 2 (R),- Fi (R 


(74) 


(75) 


(76) 


V(x) = V(R). 


R 2 v i? 2J 

Substituting ( [76| into our form of the integrals p0| we find that in our choice of gauge 
we have in fact 

Xx =h, X 2 =p 3 , (77) 

i.e. the first order integrals are again of direct geometric origin and there is no other 


possibility for them if we assume their form as in (30). 


The computation in the quantum case is essentially the same. We obtain the 

(78) 


same structure of the potentials (76) and the integrals 

A'i = L :i , 


X 2 = P 3 . 

Thus the stationary Schrodinger equation separates in the polar coordinates 

x = Rcos<j), y = Rsincj), z = z (79) 


as follows 


ip(x) = exp(im<p) exp(ikz)p(R), 


h 2 p(R) = — h 


R 


(Fi(i?) - hk) p{R) + 2 V{r)p{R) 


-A( F fi R ) + hrn ) 2 p(R)-2Ep(R). 


(80) 











Three-dimensional superintegrable systems in a static electromagnetic field 


16 


Like in the previous section let us now turn our attention towards possible 
superintegrable Hamiltonians which are integrable by virtue of the integrals 
We assume first an additional integral in the first order form 


X 3 = ^ (73 h + PI Pi) + m 3 (x). 


(81) 


i=1 


As before, the compatibility of equations (15) strongly restricts the possibilities. 
Namely, 

• if (yg) 2 + ( 7 I ) 2 7 ^ 0 then the magnetic field B(x) must vanish, i.e. this case is of 
no interest to us here; 

• if 73 = 73 = 0 then we arrive at the already discussed superintegrable case of the 
constant magnetic field and vanishing electric field. 

Thus we have not found any nontrivial Hamiltonian with magnetic Held 


superintegrable at the first order with the integrals (77). 

Next we consider the same problem with the second order integral of the form 0 
By subtracting a function of the known integrals and the Hamiltonian and using the 
relation p ■ l = 0 we can a priori set to zero the constants 

CHl, 0!14, < 233 , £*36, Q!66- (82) 

We substitute these into ( p~7[ ) and consider equations (14), (15) and (16) and their 
compatibility conditions. We find after a tedious but straightforward calculation 
whose details we are not presenting here that for nonconstant functions F\ and/or 
F '2 in (76) no second order integral <0 independent of the Hamiltonian, X\ and X 2 
can be found. 

Thus we have to conclude that the system with the potentials and field 


strength (76) is not first or second order minimally superintegrable for any nonconstant 


choice of the magnetic field B[x) and the electrostatic potential V(x). The same result 


applies also to the quantum case where only the difference between equations (16) 


and (24) needs to be considered. 


5. Superintegrability for the integrable system with integrals Li,L 2 ,L 3 
Let us now turn our attention to the case when we have three first order integrals of 


motion (32). We cannot choose among them two in involution but we easily obtain a 


second order integral 

(X) 2 = (A'0 2 + (A 2 ) 2 + (X 3 ) 2 (83) 

which is in involution with all of them. Thus assuming that we have the integrals 

X\ = l 3 +m 1 (x), X 2 =li 

we have immediately a minimally superintegrable system. 


L2 = li -1- 7712(50 


The compatibility of equations (15) for the three integrals Ai,A 2 and A 3 = 
I 2 + m 3 {x) leads directly to the following 1 -parameter family of solutions for the 
magnetic field 


B(x) = g- 


(84) 
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i.e. the only possibility is a magnetic monopole of an arbitrary strength g. The vector 
potential A is always singular at least along a halfline connecting the origin to infinity. 
We can take e.g. 


A(x) 


9 

\x\(x 2 + y 2 ) 


(y(~ 


\x\),-x{z- |x|),0) 


(85) 


which satisfies the Coulomb gauge condition X7A = 0. The functions rrij(x) are 
obtained by integrating equations (15) and up to an irrelevant additive constant read 

m i(x) = g^A. ( 86 ) 


From the condition (16) we find that the electrostatic potential V(x) must be 
spherically symmetric, 


V(x) = V(|f|). 


(87) 


Thus the classical Hamiltonian system 0 with the potentials and field strengths 
defined in (841, (85), (87) is the only system which possesses the three first order 
integrals (32) and is minimally superintegrable due to the functionally independent 
integral A 3 . Explicitly, the integrals of motion in our choice of gauge (85) read 


Ad 

A 2 


h + 9 

h + g 


- -) 

x 2 +y 2 

y{\£\ - z ) 

X 2 _|_ y2 


( 88 ) 


A3 = I3 + g, 

(A) 2 = (li) 2 + ( I 2 ) 2 + {h)~ + 2 g(l 3 + 9 )\x\ ^ , 2 ‘ 

x 2 + jr 

Since all the conditions considered are the same for both the classical and quantum 
situation, we have the same structure of minimally superintegrable system also at 
the quantum level. This rotational invariance of the magnetic monopole was already 
observed in [25] , 


Next, we shall look for an additional independent integral A 4 of the form ©> 
i.e. at most second order in momenta, which would make our system maximally 
superintegrable. That means looking at the conditions (13)-(16) for the already 
determined magnetic field B (84) and restricted electrostatic potential V(|a?|) (87) 


and establishing for which choices of V(x) an additional integral exists. We assume 
that a suitable polynomial combination of the Hamiltonian and the already known 
integrals Ai, A 2 , A 3 was subtracted from X 4 . Together with the relation 

l ■ p= 0 


we can thus set to zero the following constants aij in the integral © 


Oil, 014, OL 44, «45, 046, OL 55, 056, Ot 66- 


The compatibility conditions for equations (21) lead to the following values of the 
remaining constants a ij in © 

Ol 2 = 0, Oi3 = 0, Of 22 = 0, a 2 3 = 0, 


024 = —«15, 0 2 5 = 0, 033 = 0, 

034 = — Oi6, C*35 = —0 2 6, 036 = 0 


(89) 
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leaving three yet undetermined constants 015 , 016 , 026 - Solving the conditions (14) 
for s we find 


s ^~ 9 ^ ai5 W\ +ai6 W\ 

z x 

S 2 (x) = g ( a 26 jvjj- — 


(90) 


s 3(^) = -9 <*167^7 + 


x 


CK26] 


y 


The compatibility of the conditions (15), e.g. d x (d y m) = d y {d x m). requires that either 


• 015 , 016,026 are all equal to zero and thus also s = 0, i.e. there is no additional 
independent integral, or 

• the scalar potential must satisfy 

= ( 91 ) 

i.e. the particle moves in the Coulomb potential modified by the \x\~ 2 term 
proportional to the strength of the magnetic monopole. Solving the remaining 


conditions (15) we find that the scalar part of the integral X 4 now reads 


2 Q 

m(x) = — w (oi 6 y - a 26 x - 015 z). 
\x\ 


(92) 


The condition (16) is satisfied identically after plugging-in (90) and (91). 


The three constants 015,0161026 remain arbitrary and we have three additional 
integrals of the given form which are the components of the Laplace-Runge-Lenz 
vector modified by the presence of the magnetic monopole 

Rj = £jki (Pk + A k ) Xi — QfiL , j = 1,2,3. (93) 

\ x \ 

Of course, only one of them is functionally independent of the Hamiltonian and 
the integrals X-[. X 2 , X :i . 

The same conclusions apply also to the quantum case where the analysis is 
essentially the same, taking into account appropriate symmetrization, and the 


expressions (84), (90), (91) and (92) can be taken over literally. 


The fact that the system with the magnetic field (84) and the modified Coulomb 


potential (91) is maximally superintegrable has of course been known for long 


time (see e.g. [SHED and references therein). Here we have shown that under 
the restrictions imposed on the structure and order of the integrals there is no 
other maximally superintegrable case in this class. We notice that the restrictions 
imposed are more general than the ones under which a similar result was derived 


While it may be surprising that no modification of the isotropic harmonic 
oscillator arose in our calculation, we refer the reader to mM where it was 
demonstrated that it is maximally superintegrable but with fourth order integrals, 
not at most second, as considered here. 


6. Conclusions 


In Section [ 2 ] we derived the determining equations (13)-(161 for the coefficients of 
a general second order integral of motion (10) and discussed their compatibility 
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conditions. As in the case of a purely scalar potential {8j the coefficients hj(x),rij(x) 
satisfy equations (fl3|) that do not depend on the potentials. These equations |l3|) are 


easy to solve and imply that the leading terms in the integral X lie in the enveloping 
algebra of the Euclidean Lie algebra e 3 . As opposed to the scalar case (A(x) = 0), first 
order terms in X are not excluded (even and odd terms do not commute separat ely) . 
Also, for A(x ) 7 ^ 0 the classical and quantum determining equations differ, see (24). 
We recall that for scalar particles the classical and quantum determining equations 
and hence also the integrals of motion and the superintegrable Hamiltonians differ 
only for integrals of order N > 3 [53 55 i55 ( 1571155] . 

In general, these determining equations and also their compatibility conditions 
are difficult to solve. Indeed, if the magnetic held B(x) and the potential V(x) are not 
known the equations are nonlinear. They can be used in several manners. First of all, 
for the system to be second order integrable two such integrals must exist in addition 
to the Hamiltonian and they must commute. The leading order terms of such pairs 
of commuting integrals were classified into 11 conjugacy classes [S] under Euclidean 
transformations and this classification remains the same for B(x) 7 ^ 0 though the 
nonleading terms are different. In the purely scalar case each class corresponds to 
the separation of variables in the Hamilton-Jacobi equation and in the Schrodinger 
equation, respectively. 

This correspondence no longer holds in the presence of a magnetic held. However, 
for B(x) 7 ^ 0 it is still possible to consider each class separately and this allows 
signihcant simplihcations. The case of integrals of motion with leading terms of the 
form Pfi P| was studied in [551 . 

It is actually easier to study superintegrable systems than integrable ones since the 
conditions on the potentials are more constraining. In this article we have determined 
all superintegrable systems with at least 2 hrst order integrals and at least four 
independent integrals altogether (including the Hamiltonian). We have found the 
following superintegrable systems: 

(i) A constant magnetic held and no electric held (42). This system has 4 hrst 
order integrals (43), the Hamiltonian is a polynomial in these four. The system 
is maximally superintegrable but the fifth integral X 5 (48) is a transcendental 
function of the momenta rather than a polynomial. In classical mechanics this is 
not a problem. In quantum mechanics there are difficulties with its interpretation. 
This system is well-known to be exactly solvable [55] [56] but to our knowledge 
its maximal superintegrability has not be noticed before. The Poisson algebra of 
integrals of motion is given in (541. 

(ii) The system (591 with a periodic magnetic held and zero electric one. There are 3 
hrst order integrals pi,P 2 and ? 3 +/ 3 p 3 , the Hamiltonian H and one nonpolynomial 
integral. The classical trajectories are given in terms of elliptic functions ( 66 ) and 
examples are given on Figure [T] and Figure [2] The Schodinger equation is solved 
in terms of periodic Mathieu functions ([72]). The energy spectrum is continuous. 

(iii) The magnetic monopole with the magnetic held (84) and the scalar potential (91). 
Its second order maximal superintegrability is well-known m ■ We have shown 
that it is the only second order spherically symmetric maximally superintegrable 
system in E 3 with nonvanishing magnetic held. A fourth order superintegrable 
system is also known [53155]. 


All maximally superintegrable systems presented in this article are exactly solvable, 
both in classical and quantum mechanics. To show their superintegrability it was 
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necessary to consider nonpolynomial and nonrational integrals of motion. In a 
completely different context nonpolynomial integrals arise also for the purely scalar 
potentials, see e.g. [59114(71 fill 1421 . 

In [17] the structure of the gauge-invariant integrable and superintegrable systems 
involving vector potentials was considered in two spatial dimensions. Among other 
results it was shown there that under the assumption that the integrals are at most 
second order in momenta every superintegrable system in dimension 2 has constant 
magnetic field. However, as we have seen in Section [5] in three spatial dimensions the 
second order maximal superintegrability does not imply constant magnetic field. 
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